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Abstract. We introduce a new family of Hamiltonians with a deformed Kepler- 
Coulomb potential dependent on an indexing parameter k. We show that this 
family is superintegrable for all rational k and compute the classical trajectories 
and quantum wave functions. We show that this system is related, via coupling 
constant metamorphosis, to a family of superintegrable deformations of the 
harmonic oscillator given by Tremblay, Turbiner and Winternitz. In doing so, 
we prove that all Hamiltonians with an oscillator term are related by coupling 
constant metamorphosis to systems with a Kepler-Coulomb term, both on 
Euclidean space. We also look at the effect of the transformation on the integrals 
of the motion, the classical trajectories and the wave functions and give the 
transformed integrals explicitly for the classical system. 



PACS numbers: 03.65.FD, 02.30.K,11.30.Na 

1. Introduction 

The purpose of this article is to introduce an infinite family of classical and quantum 
systems with the Hamiltonian 



ak 2 (3k 2 
r ' 4r 2 cos 2 (|0) 4r 2 sin 2 (- 2 , 



K c =p\+ P l + V k DC , H° c = -A + V k DC (2) 

where (r, </>) are polar coordinates with < <fi < p is the linear momentum and 
A is the Laplacian on 2 dimensional Euclidean space. The superscript denotes that it 
is a deformed Coulomb potential and the subscript shows the dependence on k. This 
system was given for k = 1 in [I] and k = 2 in [2J. 

We shall show that the classical system is superintegrable for all rational values 
of k in that it allows two independent integrals of motion, besides the Hamiltonian. 
Both are polynomial in the momentum, one of second order and the other a higher 
order polynomial. We show that all bounded classical trajectories of these systems are 
closed and the motion is periodic. For the quantum system, we show that Schrodinger 
equation is exactly solvable and the energy levels are essentially the same as those 
of the Coulomb system. We will also show that these systems are related to a 
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family of superintegrable deformations of the harmonic oscillator via coupling constant 
metamorphosis. 

Superintegrable systems can be classified by the degree of the highest order 
integral of the motion, excluding the Hamiltonian. Superintegrable systems of first- 
order are directly related to Lie groups of point transformations while superintegrable 
systems of second-order are characterized by separability in multiple coordinate 
systems. Both types are considered to be well understood [3j IU El [3 [8]. The 
classification of higher-order superintegrable systems remains an open problem and 
has been a subject of much recent activity P [TU1 [HI H21 H31 [HI HS1 HH [11 HZ] ■ 

Most relevant to this paper, is the recent discovery of a family of superintegrable 
deformations of the harmonic oscillator [18] indexed by a parameter k, referred to 
by others in subsequent articles as the TTW system. It can be defined both as 
a classical or quantum Hamiltonian with potential V = uj 2 p 2 + ak 2 p~ 2 sec 2 (k9) + 
f3k 2 p~ 2 esc 2 (k9). This system has ignited much recent work on its conjectured 
superintegrability for certain values of k. Specifically, for integer k, it was conjectured 
to have an independent integral of the motion of order 2k, in addition to the second- 
order integral defining separation of variables. 

The original authors proved exact-solvability for all k and superintegrability for 
k = 1,2,3,4. Later they demonstrated the periodicity of bounded trajectories in 
the classical case for rational k [T5], supporting the conjecture. In recent papers, 
the classical superintegrability was proven for rational k [20] and the quantum 
superintegrability was proven for odd k (21 j . Most recently, a constructive proof 
of the superintegrability of the quantum system for rational k was given in 22 . 

An important tool in the analysis of superintegrable systems is the coupling 
constant metamorphosis, also referred to as the Stackel transform which maps one 
Hamiltonian to another [23] [24]. This transform is particularly useful in the study 
of integrable and superintegrable systems because it carries an associated mapping of 
the integrals of the motion. In this paper, we will see that the new Hamiltonian given 
above is Stackel equivalent to the TTW system and so the integrals of the motion, 
trajectories and wave functions will be intimately linked. 

In Sections 2 and 3, we find the classical trajectories and then the quantum 
wave functions. In Section 4 we discuss the Stackel transform and prove that systems 
with oscillator terms in the potential are Stackel equivalent to systems with Kepler- 
Coulomb terms and describe the effect of the Stackel transform on the trajectories 
or wave functions and on the integrals of the motion. In Section 5, we apply these 
theorems to determine higher-order integrals of the motion for the classical system. 



2. The Classical Trajectories 

We consider the classical trajectories of the system HJ? given by ([IK]) and prove 
that the bounded trajectories are closed and the motion is periodic. We follow the 
procedure in [25] and separate the action as S = Si(r) + S2((f>) — Et. The Hamilton- 
Jacobi equation separates as 



• dr 



A = r 2 (— ±) 2 ~Qr-Er 2 (3) 
or 



,dS 2 ^ 2 ^ ak 2 13k 1 
^ 4 cos 2 (§</>) + 4sin 2 (f - 
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and for Si and S 2 , we have 



Si(r) = ^ ~^Er 2 + Qr-Adr, S 2 (<f>) = J ^ A 
The trajectories will then satisfy 
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dS 0Si . 95 0S 2 . 

BE = BE- t = Sl ' dA = JA + ^A= 62 - (5) 
It remains to integrate these equations under the condition that the motion be 
bounded, 

< r x < r < r 2 , Erf + Qr t - A = 0, i = 1, 2. (6) 
Also, we have the requirement 

< m < cos 2 (^) < u 2 < 1, ^ M 2 + (A-/3(|) 2 + a(|) 2 K-a(|) 2 =0. (7) 

These conditions give restrictions on the choices of parameters. We summarize these 
in Table [T] below. 



Table 1. Parameter restrictions for bounded trajectories 



Restriction 


Effect 


Restriction 


Effect 


D l = Q 2 + AAE > 


r.j real 


D 2 ^(A-^(f3 + a)) 2 -^>0 


«i real 


Q > 


< r 2 


A-%\0-a\ >0 


U2 > 


A > 


< ri and 


/3 > 


Ml > 




m < cos 2 (|0) < u 2 






£ < 


r-\_ < r < r 2 


a > 


u 2 < 1 



Under the restrictions given in the table above, the solutions for ([5]) are 



-IEt-Q . 
= == h sm 



(8) 



1 . ( 2A-Qr 1 . r 2A S in 2 (^)^A+^~a) ^ 
= — 02 + ,-r arcsm( ) H arcsm( == — 3 ) (9) 



2a/I r-/DT 2fc\AA V£>2 

and hence r has the period of 2 ^—e) 3 /' 2 m ^' These trajectories are also periodic in </> 
for rational k = c/d with c, c? integer. Rewriting (j^]), 

= -2^ 2+ £±j 7r - C arcco S (^:)- rf arcco S ( Msin2( ^ ) ~* + ^ ~ " } )(10) 



2 ' rv^DT 

and using the Chebyshev polynomials defined as, 



, v , , - , sin((n+ l)arccos(a;)) 
T n (x) = cos (narccos(x)) , (7 n (x) = ; (11) 



sm arccos x 



we obtain, 



Q = -T C ( 2A ® r ) + cos(C)T d 



ry/Di 



-sin(C)i7 (1 _ 1 



2Asin 2 (|0) - A+^-Qg-a) 
75^ 



(12) 



2Asin 2 (|0) -A+^(/3-a) 



ad; 



\ 



1 - 



2Asin 2 (|0)-A+f Q8-a) 
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where C = (-2\ r ApS 2 + (c + d)/2n). Under the given restrictions on the parameters, 
the implicit function for r — r(0) given by (|12|) is well defined and periodic, with 
period r = ^tt in (f>. 

We note that these trajectories and their determining equations bear a striking 
resemblance to those obtained for the TTW system [19]. In fact, the implicit function 
determining r = r((f>) is identical under a change of variables and parameters. We 
shall see in a later section why this is so. 



3. Eigenfunctions for the Quantum System 

In this section, we solve for the wave functions of the quantum system H® c given by 
Jl|[2]) and show that the system is exactly solvable. That is, its energy values can be 
calculated algebraically and the eigenfunctions can be realized as polynomials modulo 
a gauge function (26] [27] 



We assume a solution \fr — R(r)S((j)) and separate the equation \& = E 1 ^ as 

(-d 2 r - -d r + 4 - - - e] R(r) = 0. (13) 
\ r r z r J 

l-dl + ak \ + _ a] s(^) = 0. (14) 

V 4cos 2 (4f) 4sin 2 (4f) J w y > 

We look for solutions which can be written as gauge transformations of a 
polynomial, a characteristic of exact solvability. In order to obtain such a form of 
solutions, we require that a and /3 be greater that —1/4 and rewrite a — a(a— 1), f3 = 
6(6-1). 

If we take a gauge transformation with gauge G\ — r ^A e 2r ^/-E ; the transformed 
radial equation (|13|) will have polynomial solutions if we restrict to quantized values of 
the energy E = -Q 2 {2n + l + 2y/A)- 2 . If we make a gauge transformation with G2 — 
cos(|0) a sin(|0) & then the transformed angular equation (|14l) will have polynomial 
solutions if we restrict to quantized values of the parameter A = k 2 (2m + a + 6) 2 /4. 

A set of solutions for the Schrodinger equation is given by Jacobi polynomial 
multiplied by Laguerre polynomials 



* = G 1 G 2 L 2V3 (2r V / -S)P„" 3 (-cos(fc0)) (15) 



with energy 



-Q 2 

E= 5* ff . (16) 

(2(n + km) + 1 + ka + kb) 

For a given rational k = c/d the energy levels are indexed by an integer N — dn + cm 
and their degeneracy is D — [dN/c] + 1. This coincides with the degeneracy of an 
anisotropic oscillator with frequency ratio c/d. 

These wave functions are in agreement with the solutions for the k = 1 case 
previously analyzed [?]. The quantum system is indeed exactly solvable and we 
recover the requirements A > 0, E < with a slight relaxing in the restrictions on the 
parameters a, f3 which are both required to be greater that —1/4 instead of positive as 
in the classical case. Here we see again the relation between these eigenfunctions and 
those of the TTW system which differ by the same change of variables and parameters 
as the classical trajectories. In the next section we shall see why this is the case as we 
prove some theorems which are directly relevant to our systems. 
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4. Coupling Constant Metamorphosis of Hamiltonians separable in polar 
coordinates with a harmonic oscillator term in the potential 

Consider a classical Hamiltonian H = H — EU, where H includes the kinetic energy 
and part of the potential independent of the coupling constant (-E). We can then 
write the Hamilton- Jacobi equation % = E and solve for E to obtain a new Hamilton- 
Jacobi equation % = U^(ji — E) = E. For % and % the role of coupling constant 
and the energy are interchanged. The quantum version follows by direct analogy. 
Such a transform is called coupling constant metamorphosis or Stackel transform. 
A remarkable characteristic of coupling constant metamorphosis is that there is an 
associated mapping of the integrals of the motion to the new system. Such mapping 
was given first for classical integrable Hamiltonians and later extended to quantum 
systems with second-order integrals [23j El] • More recently, the quantum transform 
was extended to higher order constants of the motion [28 . 

The theorems which define the transforms are given below with their proofs. 

Theorem 1 Given a classical Hamiltonian % = % — EU, where % is independent 
of the arbitrary parameter E, with an integral of the motion £{E). If we define the 
Stackel transform ofH and £ as H. = U~ X {H — E) and £ = £{H) then £ is an integral 
of the motion for %. 

To prove this, we use an identity for Poisson brackets given by 

d F (pi a ' 7~) 

{F(Pi,qj,f(Pi,Qj)),G} = {F{pi,qj,T),G}\ T=f{vitqj) + ^-^ — | T=/(p . >g .){/(ft, qj),G} 

where pi, qj are the conjugate position and momenta and r is a parameter. With this, 
we compute, 

{H,C} = {~(H + EU-E),£\e = *h} 

= {^(U + EU-E),£}\z =il + d E^E) \^{n , H } 
= -{U,C}^(H-E)\ s=ii 

and since 'H\^ = ^ i — E, we see that {H,£} = and we have proved the theorem. 

There is an associated theorem for quantum systems though we must make a 
further assumption about the form of the integral of motion in order to get a well 
defined integral. 

Theorem 2 Given a quantum Hamiltonian H = H — EU, where H is independent of 
the arbitrary parameter E, with an integral of the motion L = X)j=o ^-n ~2jE J , where 
Ki have degree i as differential operators. If we define the Stackel transform of H and 
L as H = U^iH - E) and L = ^.| ] K N - 2 jH d , then [H, L] = 0. 

Furthermore, if H is self-adjoint and L is self or skew adjoint, depending on the 
parity of N, with respect to d/i then H will be self-adjoint and L will have the same 
parity as L with respect to the metric U d/i. 

The proof of this theorem uses the fact that the K[s do not depend on E to show 
that for all integer j, we have 

[f] 

[L, H] = [J2 K N _ 2] E\H + EU}=0 <=> [K N - 2 j,&\ + [K N - 2j+ 2, U] = 0, (17) 
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J Lfgdp = (-1) N J fLgdp J K N _ 2j fgdp = (-1) N J fK N _ 2j gdp, (18) 

where we have extended the Ki to all integers by setting = if i < or i > N. We 
then use [K N _ 2 j, H] = U~ 1 [K N _ 2 j, H] — t/ -1 [^A^j, ^7]# to compute 

& ^ = E + [^iv- 2 ,-+2, 17]) ^ =0 (19) 

where the last equality follows by (jTTJ) . Similarly, (fT8|) implies 

jY J Kn-2iH j fgUd l i = Y^i-^f J f{HW- l K N ^ 2j Ug)Ud^ (20) 
j i 

while (HTJ) gives ^ H i U- 1 K N _2jU = V, l\ \ >,H '. Hence, the equality 
J LfgUdp = (— l) w J fLglldfi holds and we have proved the theorem. 

It is important to note that the Stackel transform usually maps one system to 
another on a different ambient manifold since there is a conformal transform of the 
metric. However, for certain forms of U, the system is still on Euclidean space, though 
with a different choice of variables; this will be the case for U = p 2 . If we assume 
further that the Hamiltonian is separable in polar coordinates, we have the following 
theorems. 

Theorem 3 Given a classical Hamiltonian % in 4- dimensional phase space, separable 
in polar coordinates and with a term in the potential corresponding to an isotropic 
oscillator, i. e. of the form 

H(p, 0)=p 2 p + \p 2 - Ep 2 + h{p) + \f 2 {B) (21) 
9 9 

where fi(r) and fi{ff) are independent of E. The Stackel transform ofH is again on 
Euclidean space and given by 

n(r, + + i/!(V27) + ±h<±) (22) 

This theorem can be directly verified by taking the Stackel transform of % as given 
in Theorem 1 

U(p, d) = j (p 2 p + ^ P 2 + h{p) + l/ 2 (0) - E 

and making the change of variables r = p 2 /2, = 29. We have the following immediate 
result. 

Corollary 1 The Hamiltonian given by h22\) is separable in polar coordinates with 
an associated integral of the motion C\ = P 2 ^ + jhi^)- Furthermore, if % has an 
additional integral of the motion then so will %. 

The first assertion can be observed from the Hamiltonian and both assertions are 
results of Theorem 1 . There is also a constructive relation between the trajectories of 
the two systems. 

Theorem 4 If H(p,9) as given in H21)) has trajectories (p(i),0(i)) which satisfy 
S± = F\ (p) — t,S 2 — F 2 ip, 9) then trajectories for 'Kir, <fi) will satisfy 

t j ' F ± iV27)dE-t, 6 2 = F 2 (V2r~,^), (23) 



Si 

dE. 
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dE 



We prove this by following the same procedure as in the previous section to solve 
for the trajectories. We separate the action as S = S\(p) + S 2 {9) — Et for H and 
§ = Si(V2r) + §2(2) ~ Et for H. By construction, the Hamilton- Jacobi equations 
are identical, the only differences being that the energy for the system H is E while 
E is the energy of the system H. Thus, we have Si(p) = Si(p) and S2{9) = S2{9). 
Therefore, the trajectories for H must satisfy 

fdS x \ dSt dS 2 

Hence, if the trajectories for T-L, (p(t),9(t)), satisfy Si — Fi(p) —t, 82 = F 2 (p, 9) then 

and so the trajectories H, {r(t), <fr(t)), must satisfy 

■ J Fx{yfo)dE-t, 62=F 2 (V2r~,^). 
We have a similar result for the quantum system. 

Theorem 5 Given a quantum Hamiltonian H(p, 9) in 2 dimensions, separable in 
polar coordinates and with a term in the potential corresponding to an isotropic 
oscillator, i. e. of the form 

H(p,0) = -~d P {pd p ) - \d 2 - Ep 2 + Mp) + \f 2 {9) (25) 
p p z p A 

where fi(p) and f%{8) are independent of E. The Stackel transform of H{p, 9) is again 
on Euclidean space and given by 

H(rA) = -\d r {rd r ) - Iflg - I + y r h(V2~r) + ^(|) (26) 

This theorem can be directly verified by taking the Stackel transform of H as given 
in Theorem 2 

H{p, 9) = \ (--d p { P d p ) - \d 2 -E + f x (p) + -i/ 2 (0) 
p \ p p p 

and making a change of variables of r = p 2 /2, <f> = 29. We have the following immediate 
result. 

Corollary 2 The Hamiltonian given by V26\) is separable in polar coordinates with 
an associated integral of the motion L\ = —d 2 , + \f2( 2 )- Furthermore, if H has an 
additional integral of the motion of the form given in Theorem 2 then so will H . 

In addition, the solutions to the Schrodinger equation for H are also solutions to the 
Schrodinger equation for H since, by construction, (H(p, 9) — E) = -^{H{p, 9) — E). 

Theorem 6 If^(p,9) is a solution to H(p,9)^(p,9) = E^(p,9) then ^(V2r, f ) will 
be a solution to H(r, 0)* ( v / 27 , f ) = E^(s/2r, §). 

If we consider the case where E = — uj 2 , fi(r) = and /a(#) = ak 2 p~ 2 sec 2 (k9) + 
(3k 2 p^ 2 esc 2 (k9) then the original Hamiltonians H, H corresponds to the TTW system, 
U\ TW ,HT TW [13]. Furthermore, if we take E = Q/2 then the transformed 
Hamiltonians H? TW , H^ TW will coincide with , Hj? c so we can use the 
superintegrability of the TTW system to prove the superintegrability of the deformed 
Coulomb system. 
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5. The Higher-Order Integral of the Motion 

We can directly apply the above theorems to the TTW system to obtain an additional 
integral of the motion for c by taking the Stackel transform of the higher-order 
constant of the motion of the TTW system. The existence of constants of the motion 
for the classical TTW Hamiltonian were recently proven for all rational k [20;. We 
would like to see the mapping of these integrals of the motion under the Stackel 
transform. In particular, we want to verify that they will remain polynomial in the 
momenta. To do this, we derive an explicit expression for the integral and show that 
these are polynomials in the parameter u 2 . 

We begin by finding an explicit expression for the constants of motion. First, we 
redefine p — e R in order to express the Hamiltonian and the integral of the motion as 

H TTW = e -2n {p 2 R+UJ 2 e ,B. + Ci) (27) 



ak 2 Bk 2 

C ^^ + ^W) + ^W) (28) 

The auxiliary functions M , Af given in [20j can be written in a modified form as 




M = — -=arccos I — x I , N = ^=arccos 




A y = A cos(2k9) - ak 2 + /3k 2 , 



B,, 



2de 



-2R 



njTTW 



A x = \/C 1 sm(2k9)pe, 

B r = 2, 



c ie - 2R PR . 



Notice that, because the quantities 

A\ +A 2 y = (A - (a + (3)k 2 ) 2 - Ak^aP, 

{nl TW f-Au> 2 c x . 



B 2 



Bi 



(29) 



depend only on Hf. , C\ and parameters, we can always multiply an integral of the 
motion by a function of these and it will still Poisson commute with T-C[ TW . 

Since the functions M and N satisfy {M,Hl TW } = {Af,Ul TW } = e~ 2R ,M—Af 
will be an integral of the motion, though not polynomial in the momenta. However, 
for rational k = c/d with c, d integer, the integral can be put into a form so that it is 
polynomial in the momenta. One such integral is 



-.(sin) 



Bl 



B 2 

v 



{^Al + Al) 



d sin (4cVA {M - Af)) 
fr- S a+ d-i 



(30) 



where 5i is if i is even and 1 when i is odd. To show that this integral is polynomial 
in the momenta, we use the identities for Chebyshev polynomials (111 |) to rewrite £2 
as 



C 



(sin) _ 



Bl + Bl 



Bn 



Bl 



A,, 



A\ + A 2 



--T r _ 



B x 



IBI + B 2 y 



B 2 

v 



U d - 



--U c -x 




A x 



+ A 2 y 



(31) 
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which can be expanded as 



C 



(sin) 



1 



E ( 2m + 1 )(-l) ms r 2ro " X Bl m+1 E( 2m H" 1 )" 1 ^" 



m=0 



m=0 



V' 



\/ -\\m Ad—2m — l a 2m+l 

2m+l K ij A * A v 



£( 2m )(-i) m ^r 2m <' 



v 



This is polynomial in the momenta because of the definition of the A' s and B' s and 
because of the parity of the Chebyshev polynomials. 

As demonstrated in the orginal paper, there is also an integral of the motion 
obtained by taking cosine instead. It can be written as 

d ros (4cVA (M - AO) 



C 



(cos) 



= JB, 



- 1 ~,. 



A i) 



(33) 



(cos) 



will 



Because of the power of \f~C\ in the denominator, the degrees of £ 2 an< ^ ^ 
differ by 1, with the lowest degree being 2(c+d) — 1. In the case of k integer, we do not 
get a symmetry operator of degree 2k but instead one of degree at least 2k + 1. In this 
case, the operator of lowest degree, namely 2k + 1 will be £2°°^ f° r even & an d ^C^ 1 "^ 
for odd k. Thus, while it has been explicitly proven that the systems are classically 
superintegrable for rational values of k, there is still no proof that the integrals of the 
motion can be written as polynomials in the momenta of degree 2k for integer k. We 
conjecture that there is such an integral £3 and it can be related to the integrals given 
above by — {£3,^1}, for fi = sin, cos depending on the parity of k. We have 
verified the conjecture for k = 1,2. 

From the equations above, we see explicitly that the integrals of motion are 
polynomial in the parameter uj 2 and hence coupling constant metamorphosis will map 
them to integrals of the motion for H'^ TW which are still polynomial in the momenta. 
To determine such integrals we must replace lu 2 = —n^ TW and perform the requisite 
change of variables. Hence, the constants of the motion for Hamiltonian "H 
be 



l TW will 



£2 (r, 



4^(727,0/2) 



2 = _-HTTW, 



(34) 



Finally, is related to n[ TW , by the parameter change E = Q/2 and so the 

integrals will need to under go the same transform. The existence of these additional 
integrals of motion proves that the classical system is superintegrable. 



6. Conclusion 



In this paper, we have presented a new infinite family of superintegrable systems 
associated with a deformation of the Coulomb potential. We have shown that the 
bounded trajectories for the classical system are periodic and that the quantum 
system is exactly solveable. We have also shown that deformed Coulomb system is 
Stackcl equivalent to the TTW system and used the classical Stackel transform to show 
that the new Hamiltonian is classically superintegrable. We mention that Bertrand's 
theorem [25] [29] (valid in n dimensions) tells us that the only spherically symmetric 
potentials for which all classical bounded trajectories are closed are the harmonic 
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oscillator and the Coulomb-Kepler potential. The potentials ar^ 1 and ar 2 are also 
the only two spherically symmetric potentials which are maximally superintegrable. 
We now see that at least for n = 2 both of these systems can be deformed into infinite 
families of superintegrable systems by adding a symmetry breaking term as in (|T|) and 
that the two families are related via coupling constant metamorphosis. 

As a tool in our analysis, we have proven some general theorems about coupling 
constant metamorphosis and its action on a class of Hamiltonians, of which the TTW 
is an example. Though the Hamiltonians , H® c are certainly novel systems, 
there is a direct relation induced by coupling constant metamorphosis between not 
only the integrals of the motion but also the trajectories and wave functions of the 
two systems. These characteristics underscore the value of the Stackel transform 
as a classifying tool. A subject of immediate interest is to try to determine which 
functions U give Stackel equivalent systems on the same manifold and also if we can 
generalize these results to higher dimensions. For example, recently a 3 dimensional 
generalization of this system was shown to be classically superintegrable for rational 
k [30]. 

Finally, it remains a further subject of research to find the closed form solutions 
of the second integrals of motion in the quantum case and prove that they can be 
chosen in a form that will admit a Stackel transform. This is the case for all explicitly 
constructed examples given in jTSJ [35] and in the general odd k case given in [5T] . 
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